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Abstract 

The  flow  phenomena  in  a  serpentine  microchannel  segment  attached  to  a  porous  transport  layer  in  a  micro  proton-exchange  membrane  fuel  cell 
is  investigated.  Due  to  the  presence  of  a  porous  transport  layer,  the  fluid  flow  for  this  configuration  exhibits  different  characteristics  compared 
with  that  through  a  simple  serpentine  channel.  The  pressure  drop  and  friction  factor  variation  in  the  channel  is  examined  for  various  values  of 
Reynolds  number  and  radii  of  curvature.  Also,  the  effect  of  variation  of  permeability  is  investigated.  There  are  two  modes  of  fluid  transport  in  this 
geometry — one  through  the  serpentine  channel  and  other  via  the  porous  media.  With  increasing  permeability,  more  fluid  is  convected  through  the 
porous  transport  layer. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Flow  phenomena  in  large-scale  fuel  cells  have  been  studied 
extensively  during  recent  years.  In  a  typical  fuel  cell  configura¬ 
tion,  fuel  and  oxidant  are  passed  through  a  serpentine  microchan¬ 
nel  attached  to  a  porous  transport  layer  (PTL).  While  passing 
through  the  channel,  the  fluid  diffuses  through  the  PTL.  The 
presence  of  a  transport  layer  attached  to  the  channel  makes  the 
flow  characteristics  significantly  different  from  those  through 
a  simple  channel.  Pharoah  [1]  studied  laminar  flow  through  a 
serpentine  channel  with  rectangular  bends  that  was  mounted  on 
a  permeable  gas  transport  layer  with  isotropic  and  orthotropic 
permeabilities.  Guvelioglu  et  al.  [2]  investigated  multicompo¬ 
nent  transport  in  proton  exchange  membrane  fuel  cells  (PEM- 
FCs)  including  electrochemical  kinetics  and  water  balance  in 
the  membrane.  These  authors  also  studied  the  effect  of  channel 
geometry  on  fuel  cell  performance.  Ramousse  et  al.  [3]  inves¬ 
tigated  heat,  mass  and  charge  transfer  in  a  single  PEMFC.  The 
independent  solutions  of  heat  and  mass  transfers  in  the  cell  were 
combined  with  a  solution  of  coupled  charge  and  mass  transfers 
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in  the  electrodes.  Hontanona  et  al.  [4]  attempted  the  optimiza¬ 
tion  of  flow  fields  in  fuel  cells.  In  particular,  the  use  of  a  grooved 
channel  distributor  and  a  porous  material  distributor  in  bipolar 
plates  was  evaluated. 

Thus,  it  can  be  seen  that  the  majority  of  the  research  work  in 
this  field  has  focused  on  species  diffusion  and  electrochemistry 
in  the  porous  transport  layer.  Also,  these  studies  are  related 
to  flow  through  rectangular  serpentine  channels.  There  is  a 
growing  interest  in  using  micro-fuel  cells  for  portable  applica¬ 
tions.  In  these  micro-fuel  cells,  the  fuel  channels  are  of  micron 
cross-section  and  Reynolds  numbers  are  typically  very  low 
(0.001-0.1).  Hence,  the  flow  through  microchannels  with  ser¬ 
pentine  bends  needs  to  be  investigated.  The  present  authors  have 
studied  flow  through  serpentine  microchannels  with  designed 
roughness  in  two  dimensions  [5]  and  three  dimensions  [6].  For 
serpentine  flow  fields  in  presence  of  a  porous  layer,  the  fluid 
has  two  different  paths  for  convection  from  the  inlet  to  the  exit 
of  the  channel.  One  is  through  the  serpentine  microchannel;  the 
second  one  is  through  the  PTL.  Depending  on  the  permeability 
of  the  porous  layer,  some  fluid  may  short-circuit  through  the 
PTL  instead  of  passing  through  the  microchannel.  Thus,  the 
convection  in  micro-fuel  cell  is  characteristically  different  from 
simple  microchannel  flow,  due  to  presence  of  a  porous  transport 
layer.  To  the  best  of  authors’  knowledge,  no  work  is  reported 
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Nomenclature 

£>h  hydraulic  diameter  of  channel  cross  section  (m) 
/  friction  factor 

Hc  channel  height  (fxm) 

//p  thickness  of  porous  layer  (pan) 

K  permeability  of  porous  transport  layer  (m2) 

Lc  length  of  straight  segment  of  channel  (fxm) 

P  pressure 

Rc  radius  of  curvature  of  channel  axis  at  bend  (jam) 
Re  Reynolds  number 

T  temperature  of  fluid  (K) 

u  velocity  component  in  v-direction  (ms-1) 

v  velocity  component  in y-direction  (ms-1) 

w  velocity  component  in  ^-direction  (ms-1) 

Wc  width  of  channel  (jam) 

Greek  letters 

e  porosity  of  porous  transport  layer 

li  viscosity  of  fluid  (N  s  m-2) 

r  shear  stress  tensor 


to  understand  such  flows  in  micro-fuel  cells.  Hence,  it  is  of 
interest  to  study  convection  in  micro-PEMFCs. 

2.  Model  for  microchannel 

The  geometry  of  the  flow  channel  in  the  micro-PEMFC  under 
consideration  is  shown  in  Fig.  1 .  The  channel  has  a  square  cross- 
section  of  Wc  =  100  [xm  x  Hc  =  100  [xm.  The  length  of  the 
straight  segment  of  channel  is  Lc  =  1000  [xm.  The  radius  of 
curvature  of  the  microchannel  at  the  bend  is  varied  from  200 
to  400  mm.  The  PTE,  located  below  the  micro-channel,  has  a 
thickness  of  Hv  =  40  |xm.  Such  flow  path  configurations  are 
commonly  observed  in  fuel  cells  in  which  a  number  of  serpentine 
channel  segments  in  series  act  as  a  passage  for  the  fuel  and  the 
oxidant.  Usually,  carbon  cloth,  carbon  paper  or  Toray  paper  is 
employed  as  the  PTE  material.  The  pore  sizes  in  such  materials 
can  range  from  submicrons  to  few  microns.  For  mass  transfer 
applications,  such  small  pore  sizes  result  in  Knudsen  diffusion, 
as  discussed  by  Litster  et  al.  [7].  Nevertheless,  here  it  is  assumed 
that  the  Brinkman-Forchheimer  extension  of  the  Darcy  equation 
is  valid  for  convection  through  the  PTE. 

3.  Governing  equations 

The  governing  equations  for  the  flow  through  the  channel  are 
the  mass  conservation  and  Navier-Stokes  equations,  which  are 
given  as: 

V(pV)  =  0  (1) 

pV(V  ■  V)  = -V P  +  iiS72V  (2) 

where  p  is  the  density  of  fluid,  which  is  air  in  this  case;  V  the 
velocity;  /x  the  viscosity;  P  the  pressure.  For  the  porous  trans¬ 


port  layer,  the  mass  conservation  and  momentum  conservation 
equations  become,  respectively: 

VfocW)  =  0  (3) 

2 

V(epW)  =  -e( VP)  +  V(er)  -  "ffv  (4) 

K 

where  e  is  the  porosity  of  the  PTE;  K  the  permeability  of  the 
medium.  The  momentum  equation  (4)  for  the  porous  media  is 
written  in  terms  of  physical  velocity,  and  simply  adds  a  source 
term  to  the  Navier-Stokes  equation  which  approaches  Darcy’s 
law  when  the  source  term  is  dominant.  With  the  inclusion  of 
the  inertia  and  the  viscous  terms,  this  equation  is  referred  to 
as  the  Brinkmann-Forchheimer  equation  and  it  always  differs 
from  Darcy’s  law  in  that  the  higher-order  terms  in  the  viscous 
term  require  the  use  of  a  no-slip  boundary  condition  at  the  solid 
walls  that  enclose  the  porous  media.  While  the  appropriateness 
of  these  additional  terms  is  debatable  [8],  this  equation  is  used 
in  a  regime  that  approaches  the  limiting  case  of  Darcy’s  law. 

In  all  cases,  the  flow  is  assumed  to  be  steady  with  constant 
properties  and  the  porous  media  are  assumed  to  be  isotropic. 

4.  Boundary  conditions 

The  channel  inlet  and  outlet  boundaries  are  shown  in  Fig.  1. 
At  the  inlet,  a  constant  velocity  is  specified  corresponding  to  a 
given  Reynolds  number.  The  inlet  velocity  is  assumed  to  be  uni¬ 
form  across  the  cross-section.  At  the  outlet,  a  constant  pressure 
equal  to  atmospheric  pressure  is  specified.  A  no-slip  boundary 
condition  is  specified  at  the  remaining  boundaries.  All  the  ex¬ 
ternal  boundaries  of  the  porous  transport  block  are  assumed  to 
be  non-permeable. 

5.  Permeability  values 


By  considering  fully  developed  Poiseuille  flow  between  par¬ 
allel  plates  and  equating  the  pressure  drop  to  Darcy’s  law,  the 
upper  limit  for  permeability  is  given  by: 


K  = 


K 


P 


12 


where  //p  is  the  thickness  of  the  PTE.  In  our  case,  this  corre¬ 
sponds  to  K  =  1.3  x  10-10  m2  for  a  40  [xm  PTE. 

An  excellent  review  of  experimental  permeability  results  for 
the  fibrous  materials  that  are  typically  used  in  PEM  PTLs  has 
been  carried  out  by  Jackson  and  James  [9].  This  paper  presents 
non-dimensional  permeability  values  as  a  function  of  the  solid 
fraction  (0=1—6)  and  shows  that  PTL  permeabilities  of  the 
order  of  10-10m2  require  fibre  diameters  of  about  45  [xm  at 
a  porosity  of  0.7  (typical  value  for  a  PTL)  with  increases  as 
the  porosity  decreases  or  as  the  desired  permeability  decreases. 
Clearly,  a  typical  fibrous  PTL  is  not  appropriate  for  a  micro-fuel 
cell  where  a  direct  scale-down  suggests  a  PTL  thickness  of  30- 
40  |xm.  It  is  possible,  however,  to  approach  the  limiting  value 
of  the  permeability  by  designing  a  new  PTL  that  is  suitable  for 
micro-fuel  cells. 

An  alternative  to  a  conventional  PTL  is  an  array  of  cylindrical 
micro  wires  connecting  the  flow-field  plate  to  the  catalyst  region. 
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Fig.  1.  Schematic  of  serpentine  microchannel. 
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where  a  is  the  cylinder  radius.  The  equation  can  be  used  to 
estimate  the  radius  required  for  a  limiting  permeability  in  the 
absence  of  sidewalls.  This  suggests  that  micro  wire  radii  would 
have  to  be  in  the  range  of  12-38  pirn  at  a  porosity  of  0.4  to 
maintain  the  permeability  between  10-11  and  10“ 10  m2. 

6.  Numerical  simulation 

The  governing  equations  are  solved  using  CFD-ACE+.  The 
geometry  is  created  and  meshed  using  the  CFD-GEOM  mod¬ 
eller.  The  solution  is  obtained  using  CFD- ACE- SOLVER  and 
the  results  are  post-processed  with  CFD-VIEW.  A  combination 
of  structured  and  unstructured  grids  is  used.  A  grid  indepen¬ 
dence  study  was  carried  out  by  successively  refining  the  mesh. 
It  is  found  that  for  number  of  cells  equal  to  447,420,  any  fur¬ 
ther  refinement  does  not  give  an  appreciable  change  in  the  so¬ 
lution.  Hence,  this  mesh  size  is  used  in  all  the  computations. 


Fig.  2.  Velocity  profile  in  serpentine  channel  at  section  Y-Y  for  Re  =  0.45, 
c  =  0.4  and  K  =  10-11  m2. 


The  density  and  size  of  these  micro  wires  can  be  set  to  obtain 
optimal  fuel  cell  performance,  and  the  desired  permeability.  The 
permeability  of  a  regular  array  of  cylinders  is  given  in  [9]  as: 


Fig.  3.  Velocity  profile  before  and  after  bend  in  serpentine  channel  for  Re  =  Fig.  4.  Variation  of  pressure  drop  with  Reynolds  number  for  serpentine  channel 
0.45,  e  =  0.4  and  K  =  10-11  m2.  with  PTL  at  e  =  0.4  and  K  =  10-11  m2. 
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The  Reynolds  number  is  varied  from  0.00045  to  0.45  and  the 
permeability  from  10-10  to  10-15m2. 

7.  Results  and  discussion 

The  surface  plot  of  the  velocity  at  section  Y-Y,  which  is  at 
the  mid-plane  along  the  channel  thickness,  is  shown  in  Fig.  2. 
At  some  distance  from  the  channel  inlet,  the  velocity  profile 
becomes  fully  developed.  Fig.  3  shows  this  fully  developed  ve¬ 
locity  profile  at  cross-sections  before  and  after  the  bend.  It  is 
observed  that  the  fully  developed  velocity  profile  continues  up 
to  the  bend.  The  velocity  profile  in  the  presence  of  a  PTL  is 
slightly  different  compared  with  a  fully  developed  velocity  pro¬ 
file  observed  in  a  flow-through  microchannel.  At  the  bend,  the 
velocity  profile  changes  due  to  curvature  of  the  channel.  After 


Fig.  5.  Variation  of  friction  factor  with  Reynolds  number  for  serpentine  channel 
with  PTL  at  6  =  0.4  and  K  —  10-11  m2. 


Fig.  7.  Variation  of  friction  factor  with  Reynolds  number  for  serpentine  channel 
without  PTL. 


the  bend,  the  profile  regains  its  fully  developed  nature  in  the 
straight  segment  of  the  channel. 

The  variation  of  pressure  drop  in  the  channel  with  Reynolds 
number  is  presented  in  Fig.  4  for  channels  with  different  values 
of  radius  of  curvature  Rc.  Here,  the  Reynolds  number  is  defined 
as: 


where  D h  is  the  hydraulic  diameter  of  the  channel.  It  is  ob¬ 
served  that  the  pressure  drop  increases  linearly  with  Reynolds 
number.  For  the  same  value  of  Reynolds  number,  the  pressure 
drop  increases  with  increase  in  the  radius  of  curvature  Rc.  This  is 
because,  as  Rc  increases,  the  effective  flow  length  also  increases. 
Hence,  more  driving  pressure  is  required  for  the  flow  to  occur  at 
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Fig.  6.  Variation  of  pressure  drop  with  Reynolds  number  for  serpentine  channel  Fig.  8.  Velocity  profile  at  interface  between  serpentine  channel  and  PTL  for 
without  PTL.  Rc  =  0.314,  e  =  0.4  and  K  =  10  ^  m2. 
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Fig.  9.  Velocity  profile  at  interface  between  serpentine  channel  and  PTL  for 
Re  =  0.314,  e  =  0.4  and  K  —  10-11  m2. 


Permeability,  m2 


Fig.  11.  Variation  of  pressure  drop  with  permeability  of  PTL  for  serpentine 
channel  at  Re  =  0.045,  Rc  =  300  |xm  and  e  =  0.4. 


the  same  inlet  velocity  (i.e.,  same  Reynolds  number).  A  plot  of 
friction  factor  against  Reynolds  number  is  given  in  Fig.  5.  The 
friction  factor  is  defined  as: 


A  P  Dh 
/_  (1/2 )pV2T 


where  L  =  Lc  +  ttRc  is  the  total  flow  length  between  chan¬ 
nel  inlet  and  outlet.  The  friction  factor  curves  for  different  val¬ 
ues  of  Rc  collapse  into  single  one,  which  indicates  that  the 
friction  factor  is  independent  of  the  radius  of  curvature  for  a 
given  channel  geometry.  The  variation  in  pressure  drop  and  fric- 
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tion  factor  for  a  channel  without  PTL  is  shown  in  Figs.  6  and 
7,  respectively.  The  plots  indicate  that  the  variation  of  these 
parameters  with  Reynolds  number  is  similar  in  nature  to  the 
flow  through  the  channel  with  a  PTL.  The  pressure  drop  and 
friction  factor  values  are  found  to  be  lower  for  channels  with 
PTLs.  The  presence  of  a  transport  layer  increases  the  avail¬ 
able  area  of  flow  in  addition  to  the  channel  flow  area.  Hence, 
this  results  in  lower  pressure  drop  for  channels  with  PTLs 
(Figs.  8-10). 

The  velocity  magnitude  along  the  interface  between  the  mi- 
crochannel  and  the  PTL  for  three  different  permeability  values  of 
10~10,  10-11  and  10“ 12  m2,  respectively.  It  can  be  seen  that  the 
velocity  increases  with  increase  in  permeability,  which  means 
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Fig.  10.  Velocity  profile  at  interface  between  serpentine  channel  and  PTL  for  Fig.  12.  Velocity  magnitude  and  vector  plot  at  bend  (section  X-X  in  Fig.  1)  for 
Re  =  0.314  and  K  =  10“12  m2.  Re  =  0.224,  e  =  0.4  and  K  =  10-'°m2. 
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Fig.  13.  Velocity  vectors  within  PTL  (section  Z-Z  Fig.  1)  for  Re  =  0.045,  e  =  0.4  and  K  =  10  10  m2. 


more  fluid  flow  occurs  through  the  PTL,  as  the  permeability  is 
increased.  At  lower  permeability,  the  PTL  behaves  as  a  solid 
plate,  which  is  confirmed  by  a  permeability  value  of  10“ 12  m2. 
The  variation  of  pressure  drop  in  the  channel  with  the  perme¬ 
ability  of  PTL  is  shown  in  Fig.  1 1 .  For  a  highly  permeable  PTL, 
the  pressure  drop  is  lower  and  increases  as  the  permeability  is 
decreased.  It  is  found  that  there  are  two  limits  for  the  pressure 
drop  and  corresponding  to  two  limiting  values  of  permeability. 
The  upper  limit  for  the  pressure  drop  is  observed  for  a  very  low 
permeability  value,  which  corresponds  to  the  case  of  a  micro- 
channel  without  a  PTL  (refer  Fig.  6),  while  the  lower  limit  is 
for  a  PTL  with  infinite  permeability.  In  this  case,  the  value  of 
the  lower  limit  is  0.38  Pa,  which  corresponds  to  a  permeability 
value  of  10-7  m2. 

Fig.  12  shows  the  velocity  at  the  centre  of  the  bend  for  the  sec¬ 
tion  X-X  (refer  Fig.  1).  It  is  found  that  a  secondary  flow  pattern 
is  present  at  the  bend  and  is  due  to  the  centrifugal  force  acting  on 
the  fluid  moving  along  the  curved  channel.  The  secondary  flow  is 
negligible  at  low  Reynolds  numbers.  It  becomes  appreciable  for 
Reynolds  numbers  in  the  range  of  4  x  10-2  to  4.5  x  10-2.  This 
secondary  flow  is  expected  to  cause  additional  pressure  drop  in 
case  of  serpentine  channels  compared  with  straight  channels.  A 
secondary  flow,  similar  to  this  one,  is  also  observed  in  serpen¬ 
tine  channels  without  porous  transport  layers  [6].  In  the  present 
case,  however  due  to  the  presence  of  porous  layer,  the  secondary 
flow  vortices  extend  into  the  PTL,  as  shown  in  Fig.  12.  This 
flow  characteristic  is  quite  different  to  the  flow  in  a  stand-alone 
micro-channel. 

The  velocity  vectors  at  a  cross-section  along  the  half¬ 
thickness  of  the  PTL  (section  Z-Z  in  Fig.  1)  are  shown  in  the 
Fig.  13.  The  serpentine  micro-channel  is  over-lay ed  on  the  figure 
to  show  the  relative  position  of  the  channel  with  respect  to  the 
PTL.  It  is  observed  that  the  fluid  flows  through  the  PTL  from  the 
inlet  to  the  exit  section  of  the  channel.  Here,  the  fluid  has  two 
paths  to  travel  between  the  inlet  and  the  outlet  of  the  channel. 
First,  fluid  flows  along  the  length  of  the  channel  and,  second,  it 
passes  underneath  of  the  channel  through  the  PTL.  The  fraction 
of  fluid  passing  through  the  PTL  reduces  as  the  permeability  is 
decreased.  It  is  also  observed  that  the  magnitude  of  the  velocity 


within  the  PTL  is  maximum  close  to  the  section  that  connects 
the  inlet  and  the  outlet,  which  corresponds  to  a  maximum  pres¬ 
sure  gradient  compared  with  other  sections  along  the  length  of 
the  channel. 

8.  Conclusions 

A  study  has  been  made  of  convection  through  a  microchan- 
nel  with  a  porous  transport  layer,  which  is  an  integral  compo¬ 
nent  of  a  micro-PEMFC.  It  is  found  that  the  velocity  profile 
in  the  channel  becomes  fully  developed  at  some  distance  from 
the  channel  inlet,  which  is  modified  due  to  the  porous  layer. 
The  pressure  drop  in  the  channel  varies  linearly  with  Reynolds 
number.  At  a  particular  Reynolds  number,  the  pressure  drop  in¬ 
creases  with  increase  in  the  radius  of  curvature.  The  variation  of 
friction  factor  with  Reynolds  number  indicates  that  the  former 
is  independent  of  the  radius  of  curvature.  The  variations  of  pres¬ 
sure  drop  and  friction  factor  in  a  channel  without  porous  layer 
are  similar  in  nature.  On  the  other  hand,  the  frictional  loss  for 
the  channel  in  the  absence  of  a  PTL  is  greater  than  that  for  a 
channel  with  PTL.  Fluid  tends  to  flow  through  the  PTL,  bypass¬ 
ing  the  serpentine  channel.  The  flow  through  the  PTL  increases 
as  the  permeability  of  the  porous  layer  is  increased.  Hence, 
convection  in  micro-PEMFCs  is  modified  by  the  inclusion  of 
a  PTL. 
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